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Abstract. Based on the results of numerical modeling, it is shown that dipole-dipole 
interactions among atoms in the active medium influences strongly the character of the 
associated superradiation. The main effect is to make the nuclear subsystem behave 
chaotically. Its strength increases with the atom density, and leads to the suppression of 
distant collective correlations and superradiation. Near correlations between the atoms are 
established, causing a confinement effect: a shielding of radiation in the active medium. 
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1. Introduction and method 

Superradiation (SR) is the cooperative radiation arising in a medium that contains a population 
inversion of excited states. Originally this effect has been stated for purely quantum systems: 
i.e. two-level atoms [1|. Experiments have confirmed this prediction [2|. Later work 
established that this phenomenon also occurs in classical systems EHH, an d that the phasing 
effect — the spontaneous origin and strengthening of correlations of originally independent 
subsystems — underlies it. In the quantum case, these are correlations among phases of 
electronic states of atoms undergoing radiative transitions; while in the classical regime 
correlations among phases of oscillations and directions of the electric dipole moments of 
atoms occur. A full account of the influence on SR of the dipole-dipole interactions among 
atoms remains incomplete (see referenses ||5]|6l|7)). 

The SR theory has been developed in several directions. There exist complementary to 
each other Schrtidinger, Heisenberg and semiclassical approaches. Each of them is applicable 
to a special area of values of the system parameters. The common methodological lack 
of these approaches is that the phasing mechanism remains off screen. The mechanism 
of the transition from casual to a phased state possesses certain spatial, time and statistical 
behaviors and its nature is not fully clear. The quantum mechanical problem of SR is rather 
complicated, for example, within the Heisenberg approach it requires to solve a system of 
nonlinear operational equations. Approximations which are used to simplify this systems 
have limited and often unclear area of applicability. Classical model of superradiation (CMS), 
where atoms are substituted by the classical Lorenz oscillators and the electromagnetic field 
is described by the classical Maxwell equations, allows to answer many difficult questions, in 
particular, the phasing mechanism. Therefore classical and quantum approaches complement 
each other. Moreover, radiation produced by pure classical system such as electrons revolved 
in magnetic field, electron clouds created in wigglers, cathode-ray lamps for microwaves, etc. 
is also SR. 

Let us consider only classical systems. First, phasing leads to the ordering of phases 
of atoms. Second, according to Earnshaw's theorem [8 , 9|, a system of point dipoles cannot 
maintain a stable static equilibrium configuration. Dipole-dipole interactions cause chaotic 
behavior that disorders their phases, and hence suppresses SR. SR arrises from a competition 
between these two opposing effects. This conclusion is inferred from the theory of non- 
uniform broadening of spectral lines for lasers ifTUl [TP . Consider now a nonlinear CMS 
lfl2l 13, i.e. a system of classical, charged anharmonic oscillators. Maxwell's equations 
describe the electromagnetic field. Next, assume that there are sufficient oscillators (N » 1), 
and they occupy a small spatial region of length L such that I « L « i, where / = n 1 ^ 3 is the 
characteristic distance between atoms, and A is the wavelength of the radiation. Each charge 
has magnitude e and mass m, and is located on the ends of springs with stiffness coefficient k, 
at coordinates r a + £ a (a = 1,2, ...,N), fixed in points r a , where there are also compensating 
charges — e. The equation of motion for the oscillators then takes the form lfl3l 

imc i m f^, \ r ab I 

HereV a = d/dr a ,r a b = r a -ri 7 ,t a b = t-r a t,/c represents the retarded time, o)q = y/k/mis 
the fundamental frequency of the oscillators, and y is the nonlinearity parameter. Substituting 
the expression 



& = b[F a (t)exp(-io)t) + F* a (t)exp(iajt)], 



(2) 
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into equation ([T) — where b represents the characteristic initial amplitude of the oscillations — 
gives 

F a + i8(\F a \ 2 -l)F a =i^V a X 

bta 

In equation (01 the second derivatives of functions F a (f) which vary slowly in comparison 
with exponents exp(+itot) are omitted, and a frequency to — too + S, 5 — 3ytoob 2 /2 is chosen. 
Note, that the case of particles rotating in a magnetic field (important in a practical sense) 
corresponds 5 < 0. For a small size system eqution (0 can be rewritten as 

r , MP i2 t\j7 3n ab (n ab F h ) - F b 1 ^ r 

F a + i6(\F a \ l - \)F a = i[32j 3 2^° Zj *' ^ 

b±a r ab b 

Where n a b = r ab /r ab , f3 = e 2 /(2mu>o), and f3o = 2e 2 u> 2 } /3mc 1 . The first term on the 
right hand side of equation represents the dipole-dipole interaction of the oscillators, 
while the second term is analagous to a 'viscosity' for the radiation in the electromagnetic 
field. Following Ref. fl2l . we shall consider one-dimensional oscillators, i.e. that dipoles 
oscillate along the x axis, and consequently, that the vectors F a are parallel to it: F a - iF a , 
i = (1, 0, 0). During a given time t we have F a (t) = p a (t) exp(iip a (t)). Hence, atoms possess a 
dipole moment that is d a (t) = e£ a = ebip a cos(<x>f + (p a ). 

The average radiation intensity of the rapidly oscillating dipoles then is 

e 2 CJ 4 b 2 sr~\ 

I = 2-1 C0S (^a - fb)- (5) 

a,b 

Thus, equation © represents a system of N oscillators, distributed arbitrarily, that can be 
solved by numerical means. A similar formalism is described in Ref. Ifl2ll ; however, dipole- 
dipole interactions are neglected. 

2. Results and discussion 

The phasing effect can be described as follows. Consider a complex plane (x,y) = 
(^R(F), 3(F)) containing points that each represent the state of an individual oscillator, 
where the distance from the origin is simply the amplitude of oscillation, and the angle is the 
phase with respect to the fundamental frequency coq. Points with to > too rotate clockwise 
around the origin; points with to < too rotate anticlockwise. 

Initially, the points are placed randomly with equal probability phases on a circle of unit 
radius p — 1 . From equation (|4j, their velocities are 

n). (6) 

b 

Here p a = (%(F a ), 3(F a ), 0), v a = p a , f = Za Pa/2, and u(p) = (0, 0, -5{p 2 - 1)), 
d(p a ,pb',r a ,rb). The latter dipole-dipole interaction term is not shown in full for reasons 
of space. Note that the vector -/ is proportional to the total dipole moment of the system 
D - ebYia Pa/2; and, to(p a ) = at t — 0. Notice also that the sign of y affects the direction 
of rotation only: changing it results in a mirror inversion without any other consequences. 
Points with positive y rotate clockwise outside the unit circle, and rotate anticlockwise when 
inside, while the opposite is true when y is negative. This symmetry, therefore, is exploited 
by choosing y > 0. 



exp(ikr ab ) 

V a x F b (t) 

r a b 
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Figure 1. Time evolution of the phase distribution of oscillators. The dotted line is a circle 
with unit radius. The number of oscillators is N = 5 X 10 3 . The concentration of oscillators 
n = 10 22 irT 3 (curve 2 on figure 4). 
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Figure 2. Time dependence of the radiation intensity for N = 5 X 10 3 (all values in arbitrary 
units). 



Having established the basis for the model, we next consider how the system evolves 
when the density of atoms n is sufficiently small that dipole-dipole interactions are negligible. 
Due to the fluctuations of density distribution of the oscillators initial phases (p a (0), the initial 
value of the vector / is not precisely zero. At t = from equation © it follows that 
dD/dt = -D/tsr, where the characteristic emission time is tsr = l/(Nfio) H E] |6] . 

Consequently from equation (0, the system responds by moving in a direction opposite 
to the dipole moment D, with a collective net velocity /. The system at time ~ t$r is 
displaced a distance ~ D(0)/(Ne) (see figure QJ). The resulting displacement moves half 
of the points outside the unit circle (p > 1), and the other half inside (p > 1). Hence, points 
outside the circle will move in clockwise orbits, while those within circulate the opposite way. 
After an interval t ~ IQtsr, the net motion results in a bunching of points on the inside of the 
circle (figureOJ)), thus the atoms emit most of their stored energy in a sharp pulse of coherent 
radiation (figure[2]i. For two-level atoms, the characteristic delay time to = tsr log given in 
IH is consistent with this. The bunch subsequently develops into a spiral-shaped distribution 
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Figure 3. Time evolution of the phase distribution of oscillators in systems with a strong 
dipole-dipole interaction. The dotted line is a circle with unit radius. Figure a and b 
corresponds to concentration of oscillators n = 8 X 10 22 irT 3 (curve 4 on figure 4 ). Figure 
c corresponds to n = 1.8 X 10 23 irT 3 (curve 6 in figure 4). 




Figure 4. Intensity of radiation (arbitrary units) for systems with different oscillator 
concentrations n (10 22 irT 3 ): 0.083, 1.0, 2.3, 8.0, 12.13, 18.38, 27.86, for 1-7 respectively. 



(figure Q]:). As it does so, the dipole moment decreases to a minimum, along with the SR 
intensity. The cycle repeats, decaying rapidly (figure |2). Oscillatory behavior is typical for 
SR in classical systems of small size fl2l . In quantum systems consisting of two-level atoms, 
SR intensity oscillations are absent (TJ. 

At high density n, dipole-dipole interactions have a significant effect. Figure[3]shows the 
outcome of equation © for large n; the initial conditions are the same as described previously. 
Notice that the points on the phase plane now move in a more chaotic manner than before. 
When n is high, dipole-dipole interactions among adjacent oscillators are strong and this leads 
to incoherence. However, SR is not entirely suppressed. In spite of the chaotic behavior 
of dipole-dipole interaction, the initial total dipole moment results in bunching of points, 
and correspondently in the SR pulse (figure |3^,b). On figure [3};, where the concentration of 
oscillators was doubled, dipole-dipole interaction suppress the bunching. 

High density systems are also complicated by collective effects. Localized groups of 
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Figure 5. Radiation intensity (arbitrary units) versus time (arbitrary units) for classical systems 
with different oscillator concentrations n (10 22 «r 3 ): 0.083, 1.0, 2.3, 8.0, 12.13, 18.38, 27.86, 
for 1-7 respectively. Case 1 is compared with the purely quantum result which varies as 
sech 2 (f - to). 



resonant atoms induce antiphase dipole moments among their neighbors. This preserves 
coherence while screening SR Q. 

The SR delay to and peak intensity I max also depend on n: increasing n makes to longer, 
and I max smaller (see figures|4]|5] and|6]l. This is a consequence of the effect of coherence on 
the collective interactions among the dipoles, which becomes weaker with increasing n. 

Unlike classical systems, quantum systems do not behave chaotically. The intensity 
varies smoothly with time as described by the following formula (TJ. 



where fi represents the form-factor of the oscillators' mutual position, and Tjv = 1 //3o is the 
characteristic emission time. This curve is plotted in figure [5] to illustrate the difference 
between the classical and quantum cases. When N is large, at t — to, equation [7] suggests 
Imax ~ N 2 - However, the CMS predicts that the exponent a = lg(I max )/ lg(N) rises to a 
peak value that is less than two, then declines as N increases (see figure |7). Experimental 
observations of SR in semiconductors exhibit similar behavior fBl . 

The results are consistent with Ref. [7]. Localized, dynamic metastable states are formed 
when the atom density n is sufficiently large. Each oscillator perturbs the motion of its nearest 
neighbors such that their relative phase differs by n. Hence, in effect each oscillator appears 
to be screened in a manner analogous to Debye shielding. This leads to confinement of 
electromagnetic fields in the active medium. 

3. Conclusions 

This study examines the phenomenon of superradiation for systems of classical nonlinear 
charged oscillators. The results of our numerical simulations show that after a characteristic 




(7) 
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Figure 6. Dependence of a maximum of radiation intensity (arbitrary units) on oscillator 
density n (in 10 22 wi~ 3 ). 
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Figure 7. Dependencies on the number of oscillators N, of (a) the ratio log 10 (/ m(1J )/ log 10 (A0; 
and (b) the peak radiation intensity. 
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delay time f , a peak in radiated power occurs, which subsequently decays in a chaotic, 
oscillatory manner, superimposed on a sech 2 (f - to) background. SR is also suppressed 
progressively with increasing oscillator density n. This behavior is ultimately a consequence 
of collective dipole-dipole interactions. These both induce incoherence among the oscillators, 
and cause a screening effect. 

Within localized regions, the individual dipoles possess correlated moments. Dipoles 
separated by sufficient large distances are nearly uncorrected. As n increases, the system 
breaks up into more of these regions. Each region emits SR impulses independently, resulting 
in the chaotic decay described above. 
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